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ABSTRACT 

We present photometric reverberation mapping of the narrow-line Seyfert 1 galaxy ESO399-IG20 performed with the robotic 15 cm 
telescope VYSOS-6 at the Cerro Armazones Observatory. Through the combination of broad- and narrow-band filters we determine 
the size of the broad-line emitting region (BLR) by measuring the time delay between the variability of the continuum and the Ha 
emission line. We use the flux variation gradient method to separate the host galaxy contribution from that of the active galactic 
nucleus (AGN), and to calculate the 5 100 A luminosity Lagn of the AGN. Both measurements permit us to derive the position of 
ESO399-IG20 in the BLR size - AGN luminosity Rblr 00 ^agn diagram. We infer the basic geometry of the BLR through modelling 
of the light curves. The pronounced sharp variability patterns in both the continuum and the emission line light curves allow us to 
reject a spherical BLR geometry. The light curves are best fitted by a disk-like BLR seen nearly face-on with an inclination angle of 
6° +3° and with an extension from 16 to 20 light days. 

Key words, galaxies: active -galaxies: Seyfert -quasars: emission lines -galaxies: distances and redshifts -galaxies: individual: 
ESO399-IG20 



1. Introduction 

The physical properties of the broad line region (BLR) in active 
galactic nuclei (AGN) have been extensively studied during the 
last three decades. Consensus has been reached that the luminos- 
ity variations of the continuum emitting region of the hot accre- 
tion disk produce variations of the broad emission lines with a 
delay due to the light travel time across the BLR. To date, the 
only method independent on spatial resolution is reverberation 
mapping ( [Blandford & McKee||1982[ |Peterson||1993| ), in which 
one measures the time delay or "echo" between the continuum 
variability and the variability observed in the broad emission 
lines. This method has been used successfully for determining 
the BLR size, kinematics and black hole mass (Mbh) of several 
Type-1 AGNs. RM provides basic Mbh estimates through the 
virial product M BH = f • R BLR • cr^/G, where G is the gravita- 
tional constant, Rblr = c • r is the BLR size, cry is the emission- 
line velocity dispersion of the BLR gas and the factor / de- 
pends on the - so far unknown - geometry and kinematics of the 
BLR ([Peterson et al. 2004 and references therein). Furthermore, 
as a widely used extrapolation of locally obtained RM results 
to higher redshift data, M BH has be estimat ed using the BLR 
Size-Luminosity 7?blr k ^agn relationship ( Koratkar & Gaskell 



1991 ;|Kaspi et al.p99g||Kaspi et al. 



Mc Gill et al.|2008||Vestergaard et al. 



20001 fWanrlel et al.|l999| 
20 lib. However, a reliable 



M B h estimation from the luminosity requires a considerable re- 
duction of the R - L dispersion present. 

As the size of the BLR ranges from a few to several hundred 
light days, it is necessary to perform observational monitoring 
campaigns ranging from months to years in order to sufficiently 
sample the time domain of the echo. Most observational RM 



campaigns have used spectroscopic monitoring of the sources, 
which however requires large amounts of telescope time. 

Recently, |Haas et aL] ( |20 1 1 \ proposed photometric reverber- 
ation mapping (PRM) as an efficient method to determine the 
BLR size through the use of broad band filters to trace the AGN 
continuum variations and narrow-band filters to catch the BLR 
emission line response. Because the narrow-band collects both 
the emission line flux and the underlying continuum, the chal- 
lenge is to extract the pure emission line light curve. In the ded- 
icated case study of 3C120, Pozo Nu nez et al. ( 2012 ) demon- 
strat ed that PRM reac hes an accuracy similar to spectroscopic 
RM flGrier et al.|2012j ). 

Another efficient appr oach to d etermine the BLR size has 
been proposed by Chelouche & Daniel (2012), through the anal- 
ysis of the difference between the cross-correlation (CCF) and 
the auto-correlation (ACF) functions for suitably chosen broad 
band filters, one filter covering a bright emission line and one 
covering only the continuum. While this method does not place 
a strict requirement on the object's redshift such that the emis- 
sio n line falls into th e narrow-band filter (as in the case proposed 
by |Haas et al.|2011|), the use of broad band filters only (as pro- 
posed by |Chelouche & Daniel |2012| ) limits the application of 
PRM to cases with a sufficiently strong emission line contribu- 
tion in the respective filter used. Despite of this handicap, this 
method has been successfully applied to determi ne the BL R size 
for one low-luminosity AGN NGC439 5 (|Edri et al.|[20T2|) and 
one high-z luminous MACHO quasar (Chelou che et al.||2012| ) 
consistent with previous spectroscopic RM results. 

A fundamental issue in both these methods is to separate the 
host-galaxy contribution from the total luminosity to calculate 
the AGN luminosity. An incorrect determination of the nuclear 
AGN luminosity (due to the contamination of the host compo- 
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nent) causes an overestimation of the linear regression slope 
a. This has been demonstrated by |Bentz et al.| ( [2009a ), who 
presented the most recent compilation of host- subtracted AGN 
luminosities for several reverberation-mapped Seyfert 1 galax- 
ies obtained through host-galaxy modeling of high-resolution 
Hubble Space Telescope (HST) images. Compared to the pre- 
vious slope of a = 0.7 they determined an improved slope of 
a = 0.5 19^o close to a = 0.5 expected from photo-ionization 



models of the BLR (Da vidson & Netzer||1979| ). A different ap- 
proach to estimate the host galaxy contribu tion is the flux varia - 
tion gradient method (F VG, |ChoToniewski| 1 98T] |Winkler|T997] ) . 
An advantage is that it does not require high spatial resolu- 
tion imaging and can be applied directly to the monitoring data. 
Recently, FVG has been tested on PRM data ( |Haas et al.|2011| 
|Pozo Nunez et al.|2012| ). These tests sho w that by using a w ell 
defined range for the host galaxy slope (Sa kata et al.||2010] ) it 
is possible to separate the AGN contribution at the time of the 
monitoring campaigns. 

ESO399-IG20 has been classified as a Narrow-Line Seyfert 
1 (NLS1) galaxy ( |Veron-Cetty & Veron|20T0| ). It has often been 
speculated that the lack of broad emission lines of - at least 
- part of the NLS1 population can be explained by face-on 
disk-like BLR geometrie s. Single-epoch spectroscopic investi- 
gations by Dietrich et al. (2005) show high-ionized gas orbiting 
at high velocities, producing strong broad profiles (FWHM(H/3) 
= 2425+121 km/s, see their Table 3 and Figures 2-5) like in 
Broad-Line Seyfert 1 (BLS1) galaxies. Furthermore, they find 
a strong host-galaxy contribution with respect to the total ob- 
served 5100A continuum flux (~ 50%). Here, we present the 
first measurement of the BLR size and the host- subtracted AGN 
5100A luminosity for ESO399-IG20. Both results allow us to 
infer its position in the BLR size-Luminos ity Rblr 00 ^gn 
gram. Additionally, pure broad-band PRM ( Chelouche & Daniel 
2012) is applied a nd compared w ith the broad- and narrow-band 



PRM technique ( [Haas et al.|2011| ). Finally, the geometry of the 
BLR, whether spherical or an inclined disk, is inferred via light 
curve modeling. 
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Fig. 1. Effective transmission of the VYSOS6 filters convolved 
with the quantum efficiency of the ALTA U16M CCD camera 
(top). The position of the narrow-band [SII] 6721 A filter (red 
line) and the simulated redshifted gaussi an- shaped Hor emis- 
sion line using the parameters obtained by Dietrich et al.| ( |200"5] ) 
(black line) are shown in the bottom panel. 



2. Observations and data reduction 

Broad-band Johnson B (4330 A), Sloan-band r (6230 A) and 
the redshifted Ha (SII 6721 ± 30 A at z=0.0249) images were 
obtained with the robotic 15 cm VYSOS6 telescope of the 
Universitatssternwarte Bochum, located near Cerro Armazones, 
future location of the ESO Extreme Large Telescope (ELT) in 
Chile. Monitoring occurred between May 5 and November 18 
of 201 1, with a median sampling of 3 days. 

ESO399-IG20 lies at redshift z = 0.0249, therefore the Ha 
emission line falls into the SII 6721 ± 30 A narrow-band filter. 
Figure [T] shows the position of the narrow-band with respect to 
the Ha emission line together with the effective transmission of 
the other filters used. The Ha line is broader than the SII fil- 
ter, however, simulations with different asymmetric and perfect 
symmetric cutting of the high- velocity line wings show negligi- 
ble systematic effect introduced in the time delay (-2%), even 
considering the fact of only using one part of the wing profile. 
This effect will be discussed in depth in a forthcoming contribu- 
tion (Bruckmann et al. 2012., in preparation). 



Data reduction was performed using IRAlf] in combination 
with SCAMP ( |Bertin|2006| ) and SWARP (Berti n etal.|2002|)rou - 
tines, in the same manner as described by |Haas et al.| ( [2012] ). 
Light curves were extracted using different apertures (5", 775, 
15" and 25") in order to compare and trace the host and AGN 
contribution. The light curves for the different apertures are 
shown in Figure [2] Additionally, for each aperture we deter- 
mined the level of variability using the fractional flux variation 
introduced by |Rodriguez-Pascual et al.| ( [l997| ): 



Vtr 2 - A 2 
if) 



(D 



where <x 2 is the flux variance of the observations, A 2 is the mean 
square uncertainty, and (/) is the mean observed flux. The vari- 
ability statistics is listed in Table [2] Selecting a small aperture 
(5") we are only considering a small portion of the total flux, 
which remains heavily dependent on the quality of the PSF, 
whereas bigger apertures (15" and 25"), enclose more contri- 
bution of the host-galaxy and the results are more sensitive to 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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Fig. 2. 2?,r-Sloan and SH-bands light curves for different apertures (5" upper left, 7.5" lower left, 15" upper right and 25" lower 
right). The light curves have been corrected for foreground galactic extinction using the values given in Table [T] 



sky background contamination resulting in larger uncertainties. 
We found that 7" 5 around the nucleus is the right size aperture 
which maximizes the signal-to-noise ratio (S/N) and delivers the 
lowest absolute scatter for the fluxes. Thus, we use this aperture 
for the further analysis of the BLR size and AGN luminosity. 
Although the light curves in Figure [2]have not yet been corrected 
for the host galaxy contribution, one can see that the fractional 
variation (F var ) and the ratio of the maximum to minimum fluxes 
(Rmax) are higher in the B-band than in the r-band for which the 
host-galaxy present a greater contribution. A different case can 
be seen in the SII light curves, which are mostly dominated by 
the contribution of the Ha emission line resulting in a higher 
amplitude of variability. As expected, for bigger apertures the B 



and r light curves become flatter with F var and R max smaller due 
to the larger contribution from the host-galaxy. 

Non- variable reference stars located on the same field and 
with similar brightness as the AGN were used to create the 
relative light curves in normalized flux units. For the absolute 
photometry calibration we used reference stars from |Landolt| 
( 2009 ) observed on the same nights as the AGN, considering 
the atmosph eric extinction for the nearby site Paranal by |Patat| 



et al.| ( |2011| ) and the recalibrated galactic fo reground extinction 
presented by Schlafly & Finkbeiner] ( |2011 ) obtained from the 



Schlegel et al.| ( |1998| ) dust extinction maps. In Table [T] we give 
the characteristics of ESO399-IG20. A summary of the photo 
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Table 1. Characteristics of ESO399-IG20 



a (2000) (1) 6 (2000) (1) z (1) 


(Mpc) 


B-V<® 
(mag) 


Ml 
(mag) 


A (3) 

(mag) 


A V) 

(mag) 


20:06:57.7 -34:32:58.0 0.024951 


102.0 


0.96 


-19.7 


0.409 


0.309 



(1) Values from NED database, (2) | Veron-Cetty & Veron| ( |2010j ), (3) |Schlafly & Finkbeiner| ( |201l| ). 
Table 2. Light Curve Variability Statistics. 







5-band 








r-band 








SH-band 






Aperture Diameter 


<f> 


<T 


F 

1 var 


Rmax 


<f> 


cr 


F 

1 var 


Rmax 


<f> 


cr 


F 

1 var 


Rmax 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


5.0 


2.24 


0.03 


0.02 


1.52 


2.86 


0.02 


0.01 


1.25 


10.56 


0.85 


0.08 


1.49 


7.5 


3.63 


0.06 


0.03 


1.45 


4.71 


0.04 


0.01 


1.24 


15.80 


1.55 


0.10 


1.44 


15.0 


5.98 


0.12 


0.01 


1.31 


8.47 


0.09 


0.01 


1.15 


23.06 


2.33 


0.09 


1.30 


25.0 


7.35 


0.14 


0.01 


1.24 


11.3 


0.09 


0.01 


1.11 


27.51 


2.42 


0.08 


1.24 



Note. Column 1 list the aperture in arc-seconds, Columns 2-5 list the mean flux, the standard deviation, the normalized excess variance and 
the ratio of the maximum to minimum fluxes for the continuum light curve. Similarly, columns 6-13 list the variability for the r and SII bands 
respectively. The standard deviation and the mean fluxes are expressed in units of mJy. 



Table 3. Summary of the photometry results for 775 aperture. 



B 

(mag) 


r 

(mag) 


SII 

(mag) 


f Z>(D 

J D total 

(mJy) 


fr a> 

J total 

(mJy) 


/snil 

(mJy) 


15.05-15.46 


14.57-14.81 


13.21-13.58 


3.45±0.06 


4.87±0.08 


17.17±0.34 



(1) 



fB total , fr total and fSll tot ai refer to the mean of the total flux ranges during our monitoring. 



metric results and the fluxes in all bands are listed in Table [3] and 
Table |4| respectively. 



3. Results and discussion 

3.1. Light curves and BLR size 

The #-band, which is dominated by the AGN continuum, shows 
a strong flux increase (about 35%) between the beginning and 
the end of May. Afterwards, the two measurements obtained in 
June reflect a more gradual increase until a maximum is reached 
at the end of July. After this maximum, the fluxes decreases 
gradually (about 20%) until the end of September, and then the 
light curves become more constant until middle of November. 
Likewise, the r-band, which is dominated by the continuum but 
also contains a contribution from the strong Ha emission line, 
follows the sames features as the 5-band light curve, although 
with a lower amplitude. In contrast to the continuum dominated 
broad-band light curves, the narrow SH-band light curve exhibits 
a flux increase which is stretched. For instance the prominent 
maximum in the B- and r- band at the end of July occurs in the 
SH-band in August, giving a first approximation for the Ha time 
delay of 15-20 days. Note that the visual inspection of the r- and 
5-band light curves does not allow an approximation for the time 
delay. 

As already discussed in previous PRM studies, the narrow- 
band contains a contribution of the varying continuum, which 
must be removed before applying cross correlation techniques 
( [Haas et al.||2011t |Pozo Nunez et al.||2012| ). In order to deter- 
mine this contribution, we used the SII and r-band fluxes, previ- 
ously calibrated to mJy, as is shown with the flux-flux diagram in 
Figure [3] The Ha line is strong contributing, on average, about 
70% of the total flux enclose in the SH-band, while the con- 




Ha =10.01 mJy ~ 72% 
^ r = 3.36 mJy - 23% 



2 3 4 5 

r flux [ mJy ] 

Fig. 3. Flux-flux diagram for the SII and r band measured using 
7.5" aperture. Black dots denote the measurement pair of each 
night. The straight red and green lines represents the average flux 
in the SII and r band respectively. The data are as observed and 
not corrected for galactic foreground extinction. 



tinuum contribution (r-band) is about 30%. Following the usual 
practice of PRM, we construct a synthetic Ha light curve by 
subtracting a third of the r-band flux (Ha = SII - 0.3 r). The 
Ha light curve was used afterwards to estimate the time delay. 
For this purpose, we used the discrete correlation function (DCF, 
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Table 4. B, r, SII and Ha Fluxes corrected by extinction. 



JD-2,450,000 




F B 
(mJy) 




F r 

(mJy) 




Fsn 
(mJy) 






Fhu 

(mJy) 




55687 


340 


2 


944 


+ 0.058 


4 


199 


+ 0.038 


12 


318 


+ 


274 


11 


058 + 


276 


55695 


316 


3 


147 


+ 0.036 


4 


279 


+ 0.052 


12 


528 


+ 


194 


11 


244 + 


200 


55698 


332 


3 


407 


+ 0.058 


4 


396 


+ 0.127 


13 


206 


+ 


281 


11 


887 + 


308 


55701 


352 


3 


631 


+ 0.051 


4 


659 


+ 0.061 


13 


561 


+ 


239 


12 


163 + 


246 


55704 


281 


3 


765 


+ 0.054 


4 


744 


+ 0.038 


14 


287 


+ 


291 


12 


863 + 


293 


55723 


359 


3 


722 


+ 0.051 


4 


763 


+ 0.060 


14 


594 


+ 


371 


13 


165 + 


375 


55726 


305 


3 


819 


+ 0.047 


4 


805 


+ 0.066 


15 


272 


+ 


161 


13 


830 + 


174 


55729 


289 




_ 


+ - 




_ 


+ - 


15 


676 


+ 


178 


15 


676 + 


178 


55745 


234 




_ 


+ - 




_ 


+ - 


15 


886 


+ 


291 


15 


886 + 


291 


55761 


137 


4 


051 


+ 0.058 


5 


073 


+ 0.064 


15 


918 


+ 


339 


14 


396 + 


344 


55765 


148 


4 


253 


+ 0.076 


5 


223 


+ 0.061 


15 


902 


+ 


291 


14 


335 + 


297 


55768 


117 


4 


149 


+ 0.080 


5 


143 


+ 0.042 


16 


273 


+ 


258 


14 


730 + 


261 


55783 


129 


3 


960 


+ 0.094 


4 


913 


+ 0.032 


16 


983 


+ 


252 


15 


509 + 


254 


55786 


133 


3 


845 


+ 0.065 


4 


885 


+ 0.080 


17 


726 


+ 


220 


16 


260 + 


234 


55795 


133 


3 


856 


+ 0.043 


4 


913 


+ 0.033 


17 


193 


+ 


194 


15 


719 + 


196 


55801 


082 


3 


772 


+ 0.043 


4 


838 


+ 0.047 


16 


774 


+ 


258 


15 


322 + 


262 


55805 


109 


3 


722 


+ 0.047 


4 


833 


+ 0.047 


16 


741 


+ 


258 


15 


291+0 


262 


55814 


109 


3 


686 


+ 0.054 


4 


861 


+ 0.088 


16 


854 


+ 


210 


15 


395 + 


227 


55817 


117 


3 


653 


+ 0.058 


4 


780 


+ 0.033 


16 


612 


+ 


194 


15 


178 + 


196 


55820 


105 


3 


566 


+ 0.047 


4 


755 


+ 0.047 


16 


322 


+ 


307 


14 


895 + 


310 


55823 


090 


3 


516 


+ 0.083 


4 


664 


+ 0.103 


16 


176 


+ 


468 


14 


776 + 


479 


55826 


094 


3 


494 


+ 0.054 


4 


645 


+ 0.052 


16 


386 


+ 


226 


14 


992 + 


231 


55829 


094 


3 


588 


+ 0.054 


4 


701 


+ 0.052 


16 


725 


+ 


404 


15 


314 + 


407 


55832 


094 


3 


577 


+ 0.065 


4 


720 


+ 0.093 


16 


531 


+ 


339 


15 


115 + 


351 


55837 


094 


3 


537 


+ 0.047 


4 


744 


+ 0.080 


16 


354 


+ 


178 


14 


930 + 


195 


55840 


102 


3 


458 


+ 0.069 


4 


683 


+ 0.042 


16 


144 


+ 


178 


14 


739 + 


182 


55843 


102 


3 


519 


+ 0.054 


4 


720 


+ 0.052 


16 


354 


+ 


339 


14 


938 + 


342 


55850 


020 


3 


588 


+ 0.058 


4 


645 


+ 0.023 


16 


192 


+ 


420 


14 


798 + 


420 


55855 


000 


3 


537 


+ 0.061 


4 


631 


+ 0.052 


15 


886 


+ 


242 


14 


496 + 


247 


55859 


059 


3 


469 


+ 0.072 


4 


598 


+ 0.099 


15 


708 


+ 


274 


14 


328 + 


291 


55869 


008 


3 


490 


+ 0.065 


4 


556 


+ 0.066 


16 


015 


+ 


436 


14 


648 + 


440 


55872 


008 


3 


519 


+ 0.051 


4 


518 


+ 0.019 


15 


983 


+ 


387 


14 


627 + 


387 


55875 


027 


3 


472 


+ 0.054 


4 


462 


+ 0.061 


16 


015 


+ 


210 


14 


676 + 


218 


55881 


016 


3 


436 


+ 0.054 


4 


556 


+ 0.061 


15 


999 


+ 


307 


14 


632 + 


313 


55884 


016 


3 


498 


+ 0.051 


4 


645 


+ 0.077 


16 


031 


+ 


500 


14 


637 + 


505 



Edelson & Krolik 1988) to cross-correlate the continuum and the 
synthetic Ha emission line, taking into account possible bin size 
dependency (Rodriguez-Pascual et al. 1989). The centroid in the 
cross-correlation of 5-band and Ha show a time delay of 18.1 
days, while the centroid in the cross -correlation between B-band 
and SH-band yields a time delay of 15.8 days. However, this 
is an expected result because of the more pronounced peak at 
zero lag due to the contamination by the continuum emission 
in the narrow-band filter. Both cross-correlation functions are 
shown in Figure [4] Uncertainties in the time delay were calcu- 
lated using the flux ran domization and random subset selection 
method (FR/RSS, |Peterson et aT1|1998b[ |Peterson et aT1[20Q4] ). 
From the observed light curves we create 2000 randomly se- 
lected subset light curves, each containing 63% of the original 
data points (the other fraction of points are unselected according 
to Poisson probability). The flux value of each data point was 
randomly altered consistent with its (normal-distributed) mea- 
surement error. We calculated the DCF for the 2000 pairs of 
subset light curves and the corresponding centroid (Figure [5}. 
From this cross-correlation error analysis, we measure a median 
lag of r cent = 18.7 +22- Correcting for the time dilation factor 
(1 + z = 1.0249) we obtain a rest frame lag of 18.2 ± 2.29 days. 

An alternative method to estimate the time delay, 
called Stochastic Process Estimation for AG N Reverberation 
(SPEAR), has been worked out recently by Zu et al. ( |2011| ). 



Through the modelin g of the AGN light curves as a damped 
random walk (DRW) ( |Zu et al.||20TT| |Zu et al.|2012| and refer- 
ence therein), this method appears to be consistent with previous 
cross-correlation techniques (e.g. |Grier et al.| ([2012 ) for spectro- 
scopic reverberation data). Using the SPEAR cod^j] on our light 
curves, we obtain a time delay of r spear = 18.4 ^ J q. After correc- 
tion for time dilation the rest frame lag r rest = 17.9 ±1.1 days, 
in agreement with the results from the DCF method. Figure [6] 
shows the SPEAR light curve models. 



3.2. Broad-band reverberation mapping approach 

Recently, pure broad-band photometric RM has been introduced 
as an efficient alternative to measure the BLR size in quasars 
(Chelouche & Daniel 2012). In this method, the emission line 
is also measured using a broad-band filter, and the removal of 
the underlying continuum is performed in the correlation do- 
main. The centroid of the lag is obtained by the subtraction 
of the auto-correlation function (ACF) between the continuum 
(represented by one broad band filter enclosing only continuum 
emission) from the cross-correlation function (CCF) between 
the continuum and the emission line (represented by one broad- 
band filter that contains a sufficiently strong emission line con- 

2 http://www.astronomy.ohio-state.edu/ yingzu/codes.html. 
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Fig. 4. Cross correlation of B and SII light curves (dotted line) 
and of B and Ha light curves (dashed line). The error range 
(±lcr) around the cross correlation was omitted for better view- 
ing. The red and black shaded areas marks the range used to 
calculate the centroid of the lag (vertical red and black straight 
lines). 
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Fig. 5. Results of the lag error analysis. The histogram shows 
the distribution of the centroid lag obtained by cross correlating 
2000 flux randomized and randomly selected subset light curves 
(FR/RSS method). The black area marks the 68% confidence 
range used to calculate the errors of the centroid (red line). 



tribution). This method has been applied for tw o objects; the 
low-luminosity AGN NGC4395 (Edri et al.|20l2) and the high- 
redshift (z = 1.72) lum inous MACHO quasar (13.6805.324) 
(Chelouche et al. |2012| ); in both cases a successful recovery 
of the time delay has been reported. Our ESO399-IG20 has a 
strong Ha emission line contributing to about 20% to the r-band 
flux (Fig. [3} and the light curves are well sampled. This makes 
ESO399-IG20 an ideal object to perform a further test of the 
pure broad-band PRM method. 

The B-band is mostly dominated by the continuum, and the 
r-band contains a strong Ha emission line. This allows us to cal- 
culate the line-continuum cross-correlation function CCF(r) = 
CCF Br (r)- ACFb(t) and to obtain the time delay (centroid). As 



Fig. 6. Synthetic Ha and continuum light curves. The solid red 
and blue lines show the Ha and continuum models estimated 
by SPEAR respectively. The red and blue area (enclosed by the 
dashed line) represent the expected variance about the mean light 
curve model. The original Ha light curve is computed by sub- 
tracting a scaled r curve from the SII curve and re-normalizing 
it to mean = 1 (red dots). The Ha light curve is vertically shifted 
by 0.2 with respect to the continuum light curve (blue dots) for 
clarity. 




Time lag [ days ] 

Fig. 7. Broad-band RM results. Top: The £-band ACF (dotted 
line) show a well defined centroid at zero lag (vertical blue line), 
while the B/ r-band CCF shows the same peak at zero lag, a sec- 
ondary peak at 3 days is also visible (vertical black line). In the 
same case as the 5-band ACF, the r-band ACF (red line) show a 
peak at zero lag, however, one additionally faint peak is clearly 
visible between 10-20 days, in comparison with the ACF of the 
5-band. Bottom: The cross-correlation yields a centroid lag of 
15.0 days (vertical black line). The red shaded area marks the 
range used to calculate the centroid. 
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Fig. 8. Results of the lag error analysis of CCF Br (r)-ACF B (r). 
The histogram shows the distribution of the centroid lag ob- 
tained by cross correlating 2000 flux randomized and randomly 
selected subset light curves (FR/RSS method). The black shaded 
area marks the 68% confidence range used to calculate the errors 
of the centroid (red line). 



shown in Figure [7] the B/r cross-correlation exhibits two differ- 
ent peaks, one peak at zero lag (from the auto-correlation of the 
continuum) and one peak at lag ~3 days. Furthermore, a small 
but extended enhancement can be seen in the auto-correlation 
of the r-band at about 10-20 days, in comparison with the ACF 
of the #-band. This feature can be interpreted as a clear contri- 
bution from the Ha emission line to the broad r-band filter. In 
fact, the cross-correlation shows a broad peak with a lag of 15.0 
days as defined by the centroid in Figure [7] To determine the lag 
uncertainties, we applied the FR/RSS method. Again from the 
observed light curves we created 2000 randomly selected sub- 
set light curves, each containing 63% of the original data points, 
and randomly altering the flux value of each data point consis- 
tent with its (normal-distributed) measurement error. We calcu- 
lated the DCF for the 2000 pairs of subset light curves and the 
corresponding centroid (Figure [5]). This yields a median lag r cent 
= 11.9 days. Correcting for time dilation we obtain a rest 
frame lag r rest = 17.5 ±3.1 days. This nicely agrees with the re- 
sults from narrow-band PRM (18.2 + 2.3 days for the DCF, and 
17.9 ±1.1 days for SPEAR) 



3.3. Modeling the geometry of the Broad-Line Region 

We have inferred the geometry of the BLR through the direct 
modeling of the results obtained from PRM. Following Welsh & 
Home (1991), the time delay function for a spherically symmet- 
ric geometry of the BLR is: 



A(r) = -(1 + sin • cos 6) 

c 



(2) 



where r is the radius of the spherical shell, and 6 the respec- 
tives angles in the spherical coordinates. While for a Keplerian 
ring/disk structure the time delay function is: 




4.0 6.2 6.3 10.5 12.7 14.B 17.0 19.2 21.3 23.5 25.7 27.6 30.0 ■ 

Inclination 



50 100 150 200 

Days after MJD 55682 (1 .5.201 1 ) 

Fig. 9. BLR disk/sphere-model. Simulated Ha light curves for a 
disk-like BLR geometry and different inclinations (4 < i < 30°) 
are shown in different colors for illustration. The black dotted 
line is the Ha simulated light curve for a spherical BLR model. A 
disk-like BLR model with inclination of 6° is able to reproduce 
the features of the original Ha light curve (black crosses). 



where r is the radius of the ring, i is the inclination (0 < i < 90°) 
of the axis of the disk with respect to the observer line of sight 
(0=face-on, 90=edge-on) and is the azimuthal angle between 
a point on the disk and the projection of the line of sight onto 
the disk. The observed continuum light curve (previously inter- 
polated by SPEAR) was convolved with the time-delay function 
for the respective disk and sphere geometry to simulate the ex- 
pected Ha light curve. The results of the simulation are shown 
in Figure [9] 

The Ha light curve for a spherical BLR configuration does 
not reproduce the observations, however a nearly face-on disk- 
like BLR geometry with an inclination between 4 to 10° provides 
an acceptable fit to the original data. To estimate the disk incli- 
nation, we performed the x 1 minimizing fitting procedure for a 
range of inclination angles (Fig [9]). The best fit yields a value of 
i = 6° ± 3° for a disk-like BLR model with an extension from 
16 to 20 light days. For such a nearly face-on BLR the velocity 
of orbiting broad line gas clouds will appear about a factor ten 
smaller, mimicking a NLS 1 . 

3.4. AGN luminosity and the Host-subtraction process 

To determine the AGN luminosity free of host galaxy contri- 
butions, we applied the flux variation gradient (FVG) method, 
ori ginally proposed by |C holoniewski (1981) and later modified 
by |Winkler et~aL] ( |1992| ). A detaile d description of the FVG 
method on PRM data is presented in |Pozo Nunez et aL] ( |2012| ), 
and here we give a brief outline. In this method the fluxes ob- 
tained through different filters and same apertures are plotted in 
a flux-flux diagram. The fluxes follow a linear slope represent- 
ing the AGN color, while the slope of the nuclear host galaxy 
contribution (including the contribution from the narrow line re- 
gion (NLR)) lies in a well defined range (0.4 < T^f < 0.53, 
for 87 



Sakata et al.1|20T0| 



A(r) = -(1 - sin i • cos 0) 

c 



'3 aperture and redshift z < 0.03, 
The AGN slope is determined through linear regression anal- 
ysis. The intersection of the AGN slope with the host galaxy 
(3) range yields the actual host galaxy contribution at the time of 
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Table 5. Total, host galaxy and AGN fluxes of ESO399-IG20 



fB total 


fBhost 


f B AGN 


f r total 


f r host 


f r AGN 


fAGN((l + z)5100A 


AL AAGN 5100A 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(10~ 15 ergs' 1 cmr 2 AT 1 ) 


(I0 43 ergs~ l ) 


3.45+0.11 


0.95+0.18 


2.50+0.21 


4.87±0.10 


2.51+0.20 


2.36+0.22 


2.66+0.21 


1.69+0.25 



AGN fluxes values fB AGN = fB total - fB host and fr AGN = fr total - fr host , with uncertainty range a AGN = {cr] otal + cr\ o J 




Fig. 10. Flux variation gradient diagram of ESO399-IG20 for 
7! 5 aperture. The solid lines delineate the ordinary least square 
bisector regression model yielding the range of the AGN slope. 
The dashed line s indicate the inte rpolated range of host slopes 
obtained from Sakata et al. ( |2010 ) for 11 nearby AGN. The in- 
tersection between the host galaxy and AGN slope (red area) 
gives the host galaxy flux in both bands. The dash-dotted blue 
lines represent the range of the AGN flux in both filters. 



Fig. 11. Flux variation gradient diagram of ESO399-IG20 for 
different apertures. The fluxes obtained through 7! 5 show an 
strong correlation, with a correlation coefficient of r c = 0.94. 
While for 5", 15" and 25" the correlation coefficients are equal 
to 0.81, 0.90 and 0.89 respectively. This significant decrease in 
the correlation supports once again our choice of 7! 5 as the best 
aperture. 



the monitoring campaign. Figure [T0| shows the FVG diagram 
for the B and r fluxes (corrected for galactic foreground extinc- 
tion) obtained during the same nights and through 7! 5 aperture. 
Additionally, FVGs were evaluated for different apertures, as 
shown in Figure [TT] As already noted in Section 2, a big frac- 
tion of the total flux is lost in the small 5" aperture, hence, the 
results are more sensitive to a possible underestimation of the 
real AGN and host galaxy contribution. Furthermore, one may 
expect that for larger apertures the fluxes will lie closer to the 
range for the host galaxy slope, however, it appears that the host 
galaxy is intrinsically strong and very blue, closer to the nucleus 
and in the outer parts. As for the analysis of the light curves we 
here use the results for the 7! 5 aperture. 

Averaging over the intersection area between the AGN and 
the host galaxy slopes, we obtain a mean host galaxy flux of 
(0.95 ±0.18) mJy in B and (2.51 ± 0.20) mJy in r. During our 
monitoring campaign the host galaxy subtracted AGN fluxes 
range between 2.29 and 2.71 and between 2.14 and 2.58 mJy 
in the B and r band, respectively. These fluxes are represented 
by the blue dotted lines in Figure [TOj From this range we in- 
terpolate the host-subtracted monochromatic AGN flux at rest- 
frame 5100A F 5l00 = 2.66 ± 0.22 • lO" 15 erg s~ l cm~ 2 A~ l . 
For the interpolation we assumed a power law spectral en- 
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Fig. 12. Rblr versus L, using data from 



Bentz et al 



(|2009a|) ; 



|Denney et aL] p010| ), [Doroshen ko et al.| p012| ) and |Grier et al. 
Q20121 ). As a red dot we included our results for ESO399-IG20. 
Shown is a zoomed portion containing ESO399-IG20. The solid 
dotted line show a fitted slope a = 0.519. ESO399-IG20 lies 
close to the expected slope. 
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ergy distribution (SED) (F v oc y a ) with an spectral index a = 
log(fBAGNlfrAGN)/ log(vfl/v r ), where v# and v r are the effective 
frequencies in the B and r bands, respectively. The error was de- 
termined by interpolation between the ranges of the AGN fluxes 
±cr in both filters. At the distance of 102 Mpc this yields a host- 
subtracted AGN luminosity at 5100AL A gn = (1-69 ± 0.25) x 
10 43 ergs' 1 . The total fluxes, host galaxy subtracted AGN fluxes 
and the AGN luminosity are listed in Table [5] 

3.5. The BLR size-luminosity relationship 

Estimates of the BLR size and host-galaxy subtracted AGN lu- 
minosity in the literature have been derived from several spec- 
troscopic RM campaigns and through host-galaxy modeling us- 
ing high-resolution images from HS T. In consequence, the rela- 
tionship between the H/3 BLR size and the luminosity (5100A) 
Rblr 00 L a (Kaspi et al. 2000| ) has been improv ed consider- 
ably with the most recent slope of a = 0.519^^ dBentzetaL 
2009a). Although to date this relationship has been corrobo- 
rated for 38 AGNs, still there exist objects with large uncertain- 
ties in both measurements. In order to improve the statistic, it 
is of interest to see the position for this new Seyfert 1 galaxy 
on the BLR-Luminosity relationship. In order to obtain the H/3 
BLR radius, we used the weighted mean ratio for the time lag 
r( Ha) : t(HJ3) : 1.54 : 1.00, obtained recently by Bentz et| 
aL] ( |2010| ) from the Lick AGN Monitoring Program of 1 1 low- 
luminosity AGN. Therefore, the Ha lag of 18.2 days translates 
into an H/3 lag of 11.8 days. Figure [T5| shows the position of 
ESO399-IG20 on the Rblr-Lagn diagram. The data are taken 
from |Bentz et"aL] ( |2009a| ) and here we include the most recent 
results for particular objects obtained from spectroscopic RM 
by |Denney et al.| ( |2010| ), |Doroshenko et al.| ( |2012| ) and |Grier et 
al. ( 20 1 2| ) respectively. 



4. Summary and conclusions 

We presented new photometric reverberation mapping results for 
the Seyfert 1 galaxy ESO399-IG20. We determined the broad 
line region size, the basic geometry of the BLR and the host- 
subtracted AGN luminosity. The results are: 

1. The cross-correlation of the Ha emission line measured in 
a narrow-band filter with the optical continuum light curve 
yields a rest- frame time delay r rest = 18.2 ± 2.29 days. 
We explored the SPEAR method, and - given that Ha con- 
tributes about 15-20% to the r-band - also the capabilities 
of pure broad-band photometric reverberation mapping. The 
SPEAR method yields a rest- frame time delay of r rest =17.9 
±1.1 days, while the pure broad-band PRM method yields 
Trest = 17.5 ± 3.1 days. The results indicate that, within the 
errors, the three methods are in good agreement. 

2. We constrained the basic geometry of the BLR by compar- 
ing simulated light curves, using the size determined for the 
BLR, to the observed Ha light curve. The pronounced and 
sharp variability features in both the continuum and emis- 
sion line light curves allow us to exclude a spherical BLR 
geometry. We found that the BLR has a disk-like shape with 
an inclination of i = 6 + 3° and an extension from 16 to 20 
light days. This nearly face-on BLR can explain the appear- 
ance of ESO399-IG20 as a narrow-line Seyfert- 1 galaxy. 

3. We successfully separated the host-galaxy contribution from 
the total flux through the flux variation gradient method 
(FVG). The average host-galaxy subtracted AGN luminos- 
ity of ESO399-IG20 at the time of our monitoring campaign 



isL AG N = (1.69±0.25)xl0 43 ^r^- 1 .IntheBLRsize-AGN 
luminosity diagram ESO399-IG20 lies close to the best fit of 
the relation. 

These results document the efficiency and accuracy of photomet- 
ric reverberation mapping for determining the AGN luminosity, 
the BLR size and the potential to constrain even the BLR geom- 
etry. 
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